Electromagnetic form factor of pion in the field theory inspired approach. 

N. N. Achaso\0 

Laboratory of Theoretical Physics, S. L. Sobolev Institute for Mathematics, 630090, Novosibirsk, Russian Federation 

A. A. Kozhevniko\Q 

Laboratory of Theoretical Physics, S. L. Sobolev Institute for Mathematics, 
and Novosibirsk State University, 630090, Novosibirsk, Russian Federation 
(Dated: January 20, 2013) 

A new expression for the pion form factor F v is proposed. It takes into account the pseudoscalar 
meson loops and the mixing of p(770) with heavier p(1450) and p(1700) resonances. The expression 
has correct analytical properties and can be used in both timelike and spacelike kinematical regions. 
The comparison is made with the existing experimental data on F v collected with the detectors 
SND, CMD-2, KLOE, and the BaBaR restricted to energies below 1 GeV. A good description of all 
four data sets is obtained. In the spacelike region, upon substituting the resonance parameters found 
in the timelike one, one obtains F w in agreement with the measurements of NA7 Collaboration. 
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I. INTRODUCTION 



The pion form factor F v is an important characteristics 
of the low energy phenomena in particle physics related 
with the hadronic properties of the electromagnetic cur- 
rent in the theoretical scheme of the vector dominance 
model [H-ll| ■ There are a number of expressions for this 
quantity used in the analysis of experimental data. The 
simplest approximate vector dominance model expression 
based on the effective 7 — p coupling cx p^A^ [3[, 
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(for notations see Sec. IIII[) does not possess the correct 
analytical properties upon the continuation to the un- 
physical region < s < Am^ and further to the spacelike 
region s < 0, nor does it takes into account the mixing 
of the isovector p-like resonances. Since, phenomenolog- 
ically, [5| gpnn/gp is not equal to unity— to be precise, 
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—the correct normalization ^(0) = 1 is satisfied by 
Eq. (jl.ip only approximately. Hereafter, a = 1/137 
stands for the fine structure constant. The formula of 
Gounaris and Sakurai @ respects the above normaliza- 
tion condition and has the correct properties under ana- 
lytical continuation. However, being based on some sort 
of effective radius approximation for the single p(770) res- 
onance, it is not suited for taking into account the mixing 
of p(770) with heavier isovector mesons. The expression 
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analogous to Eq. (11.11) based on the gauge invariant 7 — p 
coupling oc PpuFp U , 



F«{s) = 1 
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respect the correct normalization, but does not possesses 
correct analytical properties and breaks unitarity. The 
earlier expression [7J, l8| for F v takes into account the 
strong isovector mixing, but has the shortcoming that 
the above normalization condition is satisfied only ap- 
proximately, within the accuracy 20%. 

The applications of the Lagrangian of Kroll, Lee, and 
Zumino \3] to the calculations of F^ with the meson loop 
contributions in the field-theoretic context are given, in 
particular, in Refs. [§-[ll[. In particular, Ref. [10| con- 
tains the comparison of the theoretical F n with the exper- 
imental data in the spacelike kinematical region. How- 
ever, the authors of Ref. [l(| refrained from the appli- 
cation of their expression in the timelike region despite 
the fact that the high statistics experimental data col- 
lected with the detectors SND Q3] and CMD-2 [ij] were 
available at that time. 

The purpose of the present work is to obtain the ex- 
pression for the pion form factor which possesses the cor- 
rect analytical properties in the entire kinematic domain 
and takes into account the mixing of p(770) with the 
heavier resonances p(1450) and p(1700). By restricting 
the consideration to the inclusion of the pseudoscalar me- 
son loops 7r + 7r~ and KK , which admits the analytical 
treatment and is valid at energies below 1 GeV, the new 
expression is found and compared with the existing data 
on Fir collected with the detectors SND [2 CMD-2 QJ], 
KLOE [3, and BaBaR [H. 

Below, in Sec.flll the method is described by which the 
loop contributions to the vector— meson propagators are 
taken into account. The expression for the form factor 
Ftt(s) is given in Sec. IIIII Section [IV] is devoted to the 
analysis of available new experimental data on ^(s) [T2I - 
lla . Section IV] contains the discussion of the obtained 
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results. The conclusions are stated in Sec. I VII The Ap- 
pendix is devoted to the description of the method by 
which the resonance mixing is taken into account. 



II. THE LOOP CONTRIBUTIONS TO THE 
VECTOR MESON PROPAGATOR 
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Let us give some details necessary for the derivation of 
the expression for the pion form factor. They refer to the 
pseudoscalar loop contributions. For the sake of brevity, 
the notation 



Pl = p(770), P 2 = p(1450), p 3 = p(1700) 



(2.1) 



FIG. 1. The meson loop diagram contributing to both the 
diagonal polarization operator n piPi — resulting, in particular, 
in the finite width of the resonance— and the nondiagonal 
one n wPj , responsible for the pipj resonance mixing; P = 



is used hereafter for the isovector resonances involved in 
the consideration. 

The starting point is the effective Lagrangian describ- 
ing the SU(3) invariant interaction of the vector reso- 
nances with the pair of pseudoscalar mesons [161 Il7| . Re- 
stricted to the couplings of the isovector resonances pi, 
i = 1,2,3, with the pair of pions and kaons (P = n, K), 
this Lagrangian looks like 

C PiPP = i 9p ^pI ji [K-d p K+ - K + d il K~- 
K°d p K° + K°dK°] + 71-8^+ 
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The partial width of the decay pi — > PP calculated from 
the above effective Lagrangian, is 
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where s stands for the (virtual) mass squared of the de- 
caying resonance pi, and 



vp(s) = \ 1 
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is the velocity of the final meson in the rest frame of the 
decaying resonance. Applying the Cutkosky cutting rule 
to the diagram in Fig. [TJ one finds that the imaginary 
part of the diagonal polarization operator caused by the 
specific real intermediate state PP is related to the corre- 
sponding partial decay width according to the expression 



imn^ i ( s ) = 7ir PiP p( s ). 



(2.5) 



In the present work, the real intermediate states t: + tt , 
K + K~~ , and K°K° are taken into account, hence 
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The diagonal and nondiagonal polarization operators 
for the specific loop PP are calculated from the disper- 
sion integral. Here the version of this integral is defined 



which automatically provides the condition H PiPj (0) = 0, 
in agreement with the conservation of the vector current. 
To this end, the dispersion relation should be written for 
the quantity TL PiP .(s)/s. Then, one has 
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One can evaluate this dispersion integral in the unphys- 
ical region < s < 4m p, where lmIl PiPj = 0, and 
no pole is encountered. But, the integral is still diver- 
gent at s' — y co. The divergence can be regularized 
by taking the cutoff s' ma ^ = A 2 . The integration can 
be fulfilled with the change of the integration variable 
(T 2 = v p (s') = l-4m 2 p /s': 
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The logarithmic divergence can be removed by fixing 
Rel(my) = 0. The diagonal elements H PiPi = H PiPi (s) 
can be represented in the form 
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where the factor 1/2 in the second term is due to the 
flavor SU(3) relation g Pi KK = \g Pi -Kir [see Eq. (|2.2[) ] and 
that two isotopic KK modes contribute. 

The expressions for 11^^(5, my) are represented in the 
following form. Since the pion is the lightest hadron, the 
function 11^ (s, my ) looks as 
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The function 11^ (s, m^) looks different depending on the mass of the vector meson my. If my > 2m k, as is the 

case for V = p(1450) and p(1700), the expression is 
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If my < 2m,K, as is the case for V = p(770), the expres- sion is 
I 
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The function vp(s) (P = n,K) is given by Eq. (|2.4[) . over the entire s axis, looks like 
while 
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Note that the expressions Eqs. 
have the property that their real parts vanish at s = m\ 

ReTL WtK {my,my) = 0. 
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The notations are as follows. The quantity 
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III. THE EXPRESSION FOR THE PION FORM 
FACTOR 



The new expression for the pion form factor which au- 
tomatically respects the current conservation condition 
F ff (0) = 1 and possesses the correct analytical properties 
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(V = pi,2, 3,^-0 is introduced in such a way that eg^y, 
where e is the electric charge, is the transition am- 
plitude. As usual, the coupling constant gy is calculated 
from the electronic width 
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of the resonance V. The matrix of inverse propagators 



(3.4) 



is responsible for the p(770)-p(1450)-p(1700) mixing @,S, 
[l8l - [2Cl | , and A = detG. See the Appendix for more detail. 
The inverse propagators of the pi— resonance (i = 1, 2, 3) 
are 
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where the diagonal polarization operator n piPi can be ex- 
pressed through the functions IX^s, m|,) and Hk(s, Ttiy) 
described in Sec. [TTJ The nondiagonal polarization oper- 
ators are the following: 
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The quantity 023 is the dimensionless phenomenological 
free parameter. No such parameter is introduced in II P1P2 
and II P1P3 because it would result in a shift of the p(770) 
resonance peak position. See the Appendix and Refs. [2, 
1- 

The term cx H piUI in Eq. (|3.ip takes into account the 
p(770) —u(782) mixing. The basic quantities in this con- 
tribution are the following. The inverse propagator of 
the meson u(782) is taken in the form 

D u = m 2 - s - iyfsTu, (3.7) 

where the energy-dependent width 

r w = r w (s) = r w37r (s) + r W7r7 (s) + r WJrr (s) 

includes the dominant decay mode cj(782) — > 7r + 7r~7r° 
and the radiative ones. The tree pion decay width is 
represented in the form 
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where H^^s) is the phase space volume of the final 
7r + 7T~7r° state: 
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Here, m is the invariant mass of the tt + tt pair while m± 
refers to the 7r ± 7r° one: 



m 2 ± = ~(s + 3m 2 — to 2 ) ± 
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and q p7r = q(y/s, to, m ff ). Here and in what follows, 
m a ,m b ) = r^7={[s - (m a + TO h ) 2 ] x 

[s - (m a - m b ) 2 ]} 1 ' 2 (3.10) 

is the momentum of the particles a or & with the masses 
TO a or m&, respectively, in the rest reference frame of the 
decaying particle whose invariant mass is \fs. The cou- 
pling constant <7 wp7r is evaluated from the u> — > 7r + 7r~7r° 
decay width. The energy-dependent radiative width 
Tvp-y{s), where V = Pi,u), P = it, t], is related to the 

radiative width on the mass shell lyp = Typ^lrriy) in 
accord with the relation 



VP-i 



( S ) = r 



(°) 

VP7 „3 



9p( m y)' 



(3.11) 



and qp(s) = q(yfs, mp, 0) is the momentum of the pseu- 
doscalar meson P in the rest frame of the decaying vector 
meson V. The quantity 



n PlW = — jW PiU + h/ sT^ 



,(s)I> 7 (s) (3.12) 



is the polarization operator of the p(770) — w(782) mix- 
ing. The real part sn piW /m 2 is chosen in such a way 
that it vanishes at s = 0, and Tl' piU] is a free parameter. 
The contributions to lmH piul from the 777 intermediate 
state can be neglected in comparison with the irj one. If 
not fitted, the masses and partial widths of particles and 
resonances involved in the treatment are taken from the 
Review of Particle Physics [|| . 

Note that the isovector-isoscalar type of weak mix- 
ing is essential only for the p(770) — w(782) system be- 
cause it is enhanced due to the small mass difference of 
these resonances. As for other isovector-isoscalar mix- 
ings p(1450) - w(782) and p(1700) - w(782), there is no 
enhancement due to the mass proximity, and one can ne- 
glect n p2 3hJ in what follows. The coupling constant of the 
direct transition lu — > 7r + 7r~ is neglected, too. The rea- 
son for this is explained in the Appendix. See Eq. (|A8[> 
and the discussion around it. The quantities gn, g\2, gi3 
are, respectively, 
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See Eq. (|A5|) in the Appendix. 

When checking the form factor normalization (0) = 
1, one should have in mind that the pui mixing is negligi- 
ble at s = 0, because, at this energy squared, there is no 
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enhancement of the effect due to the proximity of and 
Trip. The same is true for other contributions violating 
G-parity conservation. Neglecting the above contribu- 
tions results in the correct normalization F„(Q) = 1, if 
one takes 
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Indeed, the mixings due to strong interactions n piPj . van- 
ish at s = 0, and F„(0) reduces to the above sum. This is 
the reason for the s in front of 023 in Eq. Q3.6p . The com- 
parison of the new expression Eq. (13.11) with the latest 
experimental data [12Hla | obtained in e + e~ annihilation 
is presented in the next section. 



y/s < 1 GeV, because, at the first stage of the study, 
the proposed expression for the polarization operator is 
restricted to include only tt + tt~ and KK loops. 

The original e+e" -> tt+tt" data of the SND, CMD-2, 
and KLOE Collaborations are presented in two distinct 
forms. The first one is the form factor with the vacuum 
polarization effect included. The BaBaR Collaboration 
does not present their results in this form. The second 
form is the so-called bare cross section. This quantity 
is undressed from the vacuum polarization effects, but 
includes the hnal state radiation. All four groups present 
their data in this form. For the purpose of uniformity of 
presentation, the analysis of the present work refers to 
the bare cross section 



IV. THE DATA ANALYSIS AND RESULTS 
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The experimental data on the reaction e + e — > 
collected by the collaborations SND jlj], CMD-2 [l3|. 
KLOE [IH, and BaBaR [ll| are chosen for the analysis 
in the framework of the field-theory-inspired approach to 
the pion form factor presented in this work. As for the 
BaBaR data set, we restrict ourselves by the points with 



where F w (s) is given by Eq. p.l[) , 

q„(s) = y/sv 1T {s)/2 

is the momentum of the final pion, and the function a(s) 
allows for the radiation of a photon by the final pions. In 
the case of the pointlike pions, it has the form [12, |22T - [25| 
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Here, v n = v 7r (s) is given by Eq. (|2.4p . and 

. . f x , ln(l - t) 

Li 2 (a:) = - / dt-±- 

Jo t 

First of all, no fit with the single p(770) resonance con- 
tribution, based on Eq. (13. ip in which both <7 P27r7r arL d 
9p 3 nir are set to zero, is capable of satisfactory descrip- 
tion of all four data sets, even with the pu> mixing effect 
being taken into account. Although the formula with the 
single resonance works well in the put resonance region, 
the curve at the far-right shoulder of the p(770) resonance 
peak does not follow the data points. 

Taking into account the resonance ^2, but with the 
neglect of the ps one, results in a rather poor fit, too. 
This is because the normalization condition -F^O) = 1 
reduces, in this case, to the rather restrictive sum rule 

.9pi7T7T 9p2TTT , 

"T ^ , 

9pi 9p2 

which fixes completely the P2 contribution to the e + e _ — > 
7r + 7r~ reaction amplitude in a way that forbids the suc- 
cessful fit. Specifically, the ratio g P2 -mt / 9p 2 turns out to be 



too small due to the fact that the universality condition 
9p\ttixI 9p x ~ 1 is satisfied for the couplings of p(770). See 
Eq. (II. 2[) . Hence, the pi resonance contribution turns out 
to be smaller than necessary for reconciling the calcula- 
tions with the data. The third resonance p$ = p(1700) is 
required in order both to preserve the approximate uni- 
versality condition and to allow a freedom in the variation 
of the P2 = p(1450) couplings. 

Free parameters, which should be determined from 
comparison with the existing data [l2|-[l5j , are the masses 
of the resonances p(770) and w(782), the coupling con- 
stants g Pl 2 3 _ i . 7rTr of the resonances pi,2.3 with the 7r + 7r~ 
state, the coupling constants g Pl 2 and parametriz- 
ing the pi,2,3 and w(782) leptonic decay widths [see 
Eq. (|3.3p ]. and the real part of the polarization oper- 
ator of the p(770) - w(782) mixing W piUJ . Note that 
g P3 is not free but should be determined from the sum 
rule Eq. (|3.13[) . At last, there is the parameter a 2 3 [see 
Eq. (|3.6p ] that defines Ren p2P3 . Since we restrict our 
analysis to the energy range below 1 GeV, the masses of 
the resonances p(1450) and p(1700) are fixed to, respec- 
tively, m P2 = 1.45 GeV and m P3 =1.7 GeV. 
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FIG. 2. The bare cross section, Eq. (|4.ip . calculated with the 
resonance parameters obtained from fitting the SND data [T3 ] 
listed in Table [I] Experimental points are from Ref. [l2| . 
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s 1 " [GeV] 

FIG. 3. The same as in Fig. [21 but evaluated with the parame- 
ters obtained from fitting the CMD-2 data pl| • Experimental 
points are from Ref. [l3T ]. 



So, the total set of free parameters is 



9p 3 iT-n, «23- 



(4.3) 



Their obtained values, found from fitting the bare cross 
section Eq. (|4.1[) . side- by-side with the corresponding x 2 
per number of degrees of freedom, are listed in Table 
U separately for the four independent measurements of 
SND Q3], CMD-2 [ij], KLOE Q, and the BaBaR data 
(llj restricted to the low-energy range y/s < 1 GeV by 
the reason explained earlier. The bare cross section eval- 
uated with the parameters of TableJJ is compared with 
the SND [3, CMD-2 [3, KLOE [lj, and BaBaR p| 
data shown in Figs. [5J El HI and 03 respectively. 

As far as the specific values of the obtained parame- 
ters in the Table Hare concerned, those corresponding to 
the /?(770) — u(782) resonance system agree satisfactorily 
for all four experiments |T2h15|. The agreement of the 
coupling constants of the resonances p(1450) and p(1700) 
is poor but, taking into account the large uncertainties 
in their determination, is not crucial. This is justifiable, 
because the energy range y/s < 1 GeV is not a proper 
place for extraction of the coupling constants of the above 
resonances. The widths of p(1450) and /o(1700), in their 
respective energy ranges, are known to be saturated by 
the complicated final states pirn, lott, etc., not the 7r + 7r~ 
one Taking into account these decay modes is nec- 
essary at energies ^fs > 1 GeV. Unfortunately, taking 
into account the real parts of the polarization operators 
arising due to the mentioned complicated states is hardly 



1000 



600 




s [GeV ] 

FIG. 4. The same as in Fig. [2] but evaluated with the param- 
eters obtained from fitting the KLOE-2010 data [HI]. Exper- 
imental points are from Ref. [lil ]. 
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TABLE I. The resonance parameters found from fitting the data from SND [T^, CMD-2 [T3], KLOE10 0, and the BaBaR 
data [IB] restricted to the energies \fs < 1 GeV. 



parameter 


OlN U 


KjNLLJ-Z 




DaDa.rv 


m pi [MeV] 


773.76 ± 0.21 


774.70 ± 0.26 


774.36 ± 0.12 


773.92 ± 0.10 




C 700 j_ n nnc 
tj. I i?0 _l_ u.uuu 


c 70c 1 n nno 
<j. / 00 _l_ u.uuo 




C 7OC J_ A QQ4 


J Hi 


5.130 ±0.004 


5.193 ±0.006 


5.242 ±0.003 


5.167 ±0.002 


m u [MeV] 


781.76 ±0.08 


782.33 ± 0.06 


782.94 ±0.11 


782.04 ±0.10 




17.13 ±0.30 


18.43 ± 0.47 


18.27 ±0.45 


17.05 ±0.29 


10 3 n^ [GeV 2 ] 


4.00 ±0.07 


3.97 ±0.10 


3.98 ±0.09 


4.00 ±0.06 




0.71 ±0.35 


0.79 ±0.26 


0.019 ±0.004 


0.21 ±0.04 


5 P 2 


8.0 ±4.4 


7.6 ±3.4 


0.22 ±0.07 


4.0 ± 1.0 


3P37T7T 


0.20± ;f 7 


0.76 ±0.75 


r)55+ ' 088 

U.UOJ_ 043 


011+ ' 479 
u.uii_ 007 


0,23 


0.002 ±0.011 


-0.016 ±0.057 


-0.014 ± 0.040 


-0.0005 ± 0.0009 


X 2 /Nd.a.f. 


54/35 


34/19 


87/65 


216/260 
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FIG. 5. The same as in Fig. [2] but evaluated with the pa- 
rameters obtained from fitting the BaBaR data [lj| restricted 
to the energies \fs < 1 GeV. Experimental points are from 
Ref. O. 



possible in closed form. In addition, the corresponding 
dispersion integrals diverge much more strongly than in 
the case of the 7r + 7r _ and KK intermediate states con- 
sidered in the present work. In the meantime, the small 
values of g P2 ^ir, in comparison with 5 Pl7r ir, obtained in 
the present work upon neglecting the pirir, unr, etc., decay 
modes at y/s < 1 GeV, agree with the earlier conclusions 
0, H| inferred from the analysis in which the above decay 
modes were included. Note also that 023 is compatible 
with zero. 



DISCUSSION. 



An important check of the expression for the pion form 
factor Eq. (|3.1[) and the consistency of the fits is the con- 
tinuation to the spacelike region t < accessible in the 
scattering processes. To this end, one should take the 
branch with s < in lT^^m^) [see Eqs. (f278|) . ([2lf]) . 
and (|2.10l) ] and replace s —> t. Having in mind that the 
p(770) — w(782) mixing in the region t < is negligi- 
bly small one can calculate F w (t) in this region. The 
results are shown in Fig. [6l where the comparison with 
the NA7 data [2lJ is presented for all four fits considered 
in the present work. We emphasize that the data [2l[ 
are not included to the fits. Hence, a good agreement, 
demonstrated in Fig. [6l makes the evidence in favor of 
the validity of Eq. (|3.1[) for the pion form factor. 

Using the resonance parameters of Table HI one can 
calculate, in particular, such important characteristics as 
the charged pion radius , defined as the square root of 
the root-mean squared radius, 



of the spherical symmetric electric charge distribution 
FAq) = J d 3 rp(r)e^ r « F ff (0) - 

^ J p( r )r 2 d 3 r = F 7r (0)+ t -(r 2 ), (5.1) 
where t = —q 2 - One gets 
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,dFjf) 



¥0- 



(5.2) 



Evaluating ?v with the parameters of TableHl one obtains 
the results presented in the first row of Table [II] For 
comparison, the averaged value of the pion charge radius 
cited by the PDG [| is rv = 0.672 ± 0.008 fm. 

If one considers the single p(770) resonance, then its 



inverse propagator near s 



in 



can be represented as 



s + (m 



8 



TABLE II. The pion charge radius, tv, Eq. 
(with the superscript phys), the bare ones 
resonance parameters of Table |U 



(|5.2[) . the renormalization constant, Z p , Eq. ()5.4[) . the "physical" partial widths 
(without the superscript), of the decay p(770) and w(782), evaluated with the 



parameter 



TV [fm] 

T pin7T {m 2 pl ) [MeV] 
rff«) [MeV] 



SND 



CMD-2 



KLOE10 



BaBaR 



- P } 7T7T 

Ep 1 ee(^ip 1 
p(phys) ^__2 

rZ e (mi) [keV] 



keV] 
,,,, e (m pi ) [keV] 



0.635 ± 0.054 
0.9273 ± 0.0003 
139.93 ± 0.29 
150.90 ± 0.31 
6.56 ±0.01 
7.07 ±0.01 
0.59 ±0.02 



0.646 ± 0.059 
0.9277 ± 0.0002 
139.54 ± 0.39 
150.42 ±0.42 
6.41 ± 0.01 
6.91 ±0.01 
0.51 ± 0.03 



0.668 ± 0.039 
0.9279 ± 0.0002 
139.12 ±0.29 
149.92 ± 0.31 
6.29 ±0.01 
6.78 ±0.01 
0.52 ± 0.03 



0.668 ± 0.053 
0.9277 ±0.0001 
139.34 ±0.19 
150.20 ±0.20 
6.47 ±0.01 
6.97 ±0.01 
0.60 ±0.02 
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FIG. 6. The pion form factor squared in the spacelike region 
evaluated using the resonance parameters of Table [I] The 
labels of the theoretical curves correspond to the columns of 
Table U The experimental data NA7 are from Ref. [2l| . 



iVsT pl ^(s). (5.3) 

The behavior of Ren piPl (s) is shown in Fig. [7] 
Comparing Eq. (|5.3p with Eq. (|l.lj) one can see that one 
should make the renormalization 



9 P i 



Z- 1/2 a 



zl /2 g P1 , 



where 



1 + 



rfRen piPl (.s) 
ds 



(5.4) 



in order to reduce Eq. (|5.3I) to the conveniently used form 
with m pi being the physical mass of the resonance. This 
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FIG. 7. The energy dependence of Rell pipi (s) for both tt + tv 
and K K loops. 



results in the renormalization of the 7r + 7r and e + e par- 
tial widths of the p(770): 



r 
r 

x pi ee 



-p(phys) _ x P17T7T 

P17T7T '/ 
Z P 



n(phys) 
piee 



r 



piee 



z„ 



(5.5) 



The numerical values of the renormalization constant Z p 
are given in Table [TTl side-by-side with the 7r + 7T - and 
e + e~ partial widths of the p(770). One can see that 
Z p brings the "bare" widths (without the superscript 
"phys") closer to the values T p7rTr = 149.1 ± 0.8 MeV and 
r pee = 7.04 ±0.06 keV cited in the Review of Particle 
Physics [5]. 

Another important characteristic of the low-energy 
hadronic physics is the phase shift 5\ of 7T7T scattering in 
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1,0 



and KK loops to both diagonal and nondiagonal polar- 
ization operators. In principle, other intermediate states 
could be taken into account, at least numerically. How- 
ever, heavier isovector resonances p(1450) and p(1700) 
are known to have other decay modes besides 7r + 7r~ and 
KK, such as ujtt, a\ir etc. The treatment should include 
the energies \fs < 2 GeV where the coupling constants 
with the above states could be determined. No data exist 
on these decay modes of the quality comparable with the 
7r + 7r~ data [12 - fl5l |. Hence, at present, the restriction to 
the domain y/s < 1 GeV and to the pseudoscalar loops 
seems justifiable. 
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Appendix A: The finite width and the resonance 
mixing 



FIG. 8. The phase shift Si of tvtt scattering. The data are, 
resp ectively, Protopopescu et al. [26| and Estabrooks et al. 
[271 ] . The curves corresponding to the parameters obtained 
from fitting the SND, CMD-2, and KLOE data are not shown 
because they coincide with the curve evaluated using the pa- 
rameters from the fit of the BaBaR data, shown here. 



the vector-isovector channel with the quantum numbers 
of p(770). At energies below the ujtt and KK produc- 
tion thresholds, 5{ is given by the phase of the pion form 
factor 



S\ = arctan 



Imi^ 
ReiV 



(5.6) 



where F v is given by Eq. (|3.1|) upon neglecting the contri- 
bution of pu mixing oc H piul . The plot of 5{, obtained us- 
ing parameters extracted from fitting the low-energy por- 
tion of the BaBaR data [l5 j. i s shown in Fig. [HJ where the 
comparison with the data |26l. [27j is presented. Note that 
the resona nce p arameters, extracted from three other sets 
of data [E2h14| . result in the curves for 8\ coincident with 
that shown in Fig. [5] Having in mind that the data on the 
phase shift were not included in the fits, the agreement of 
the calculated 5\ with the measured one is satisfactory. 



VI. CONCLUSION 



It is shown that the new formula for F„(s), Eq. (|3.ip . 
gives a good description of the latest experimental data 
[12Hl5| ] on the production of the ir + ir~ pair in e + e~ an- 
nihilation at y/s < 1 GeV. In this low-energy domain, 
one can restrict oneself by the contribution of the tt + it~ 



Some details necessary for taking into account the fi- 
nite width effects and the resonance mixing are given 
in this Appendix. The meaning of the diagonal polar- 
ization operator Hrr(s) is that it modifies the inverse 
bare propagator of the resonance R with the mass tur, 



D$\8) = D<*=m» R 
1 1 



- s, in the following way: 
1 , . 1 



Dr(s) 



D 



(0) 

R 
1 



D 



nijfl(s)— m + 



D 



(0) 

R 



1 



D 



(0) 



n jRjR (s) T ^yn jRjR (s) 



D 



D 



(0) 



1 



D$ - T1rr(s) 



'R 

In particular, this formula takes into account the finite 
width effects 

D R (s) = m\ - s - Ren RR (s) - i\/sT Rvv (s). (Al) 

In principle, the mixing of the isovector resonances 
p(770), p(1450), and p(1700) can be strong, especially 
because of the common decay modes, for example, the 
7r + 7r~ one. It can be taken into account in the field- 
theory-inspired approach based on summing to all orders 
of the loop corrections to the bare propagators of vector 
mesons 0, H, EE E3 ■ The term "bare" means that the 
propagators are not distorted by the mixing. The scheme 
can be demonstrated by taking the two-resonance mixing 
as an example [2(| ■ It reduces in this case to the following 
replacements: 



1 1 

Dr'^Dr' 



1 1 1 

7^— HaR/ — — Urri — — 

Ur Dri JJr 



D 



R' 



DrDr, 



n 2 

^RR' 



RR ■ 
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1 



1111 

L>R' MR' Ur L>r' 



D, 



D R D R , - n 2 



= (G~ 



n 



RR' 



RR' 



TIrr 



DrDr, - IP RR , 



(G- v 



The matrix 



G 



Dr -Urr, 
Urr, Dr, 



IR'R' 



RR' ■ 



is the matrix of inverse propagators in the two-resonance 
case. Let us take for a moment just this case, R = p\ and 
R' = p2, in order to clarify the effect of the mixing on 
the resonance position. Neglecting for a moment the pw 
mixing which is taken into account below, one can write 
the pion form factor as 



F-n G?7Pi ' 9*1 P2 
1 



D 

n 



P2 
P1P2 



n 



P1P2 

D P1 



9pitttt 

9 P2TZ1T 



n n — n 2 

LJ Pl 1J P2 LL pip 2 



(A2) 



In the vicinity of the pi resonance position, s — > m 2 , 



Eq. (|A2I) can be represented in the form 



m p! - 5 - n piPl(^)- 



n 2 

- m- — 



(m 2 ) 
v PI ' 



(A3) 



where, in accord with the adopted definition, 
ReII l01 p 1 (m 2 ) = 0. One can see from Eq. (|A3[) 
that there is a shift in the p\ resonance peak position 
due to the mixing of p\ with the resonance pi\ 



Am 



pi 



-Re- 



n 2 „ „ (ml " 

— "2 \ pi , 



P1P2 ' 



Re[n 2 1P2 ( 



pi 

•v,2 



npsps^pj 



(A4) 



where we neglect n p2P2 (TO pi ) in comparison with the 

Indeed, using the 



mass difference squared m 
plots in Fig. [7j the relation 



P2 



n 



9 P 2 



pi ■ 



n„ 



P2P2 \ J ix pipn 

ypi7T7r / 

Eq. Eq. pTfj) . and # P27r7r « 0.8 (see Table U, one 

obtains the estimate 

IWm|) < (0 2 + 15j)xlQ -3 

m p2 - m 'Pl 

In the case of the well-studied resonance pi = p(770), it 
is natural to expect that the visible peak position with 
a good accuracy coincides with the bare mass m pi . This 
follows from the definition Re±I plPl (m pi ) = adopted in 
the present work. In order to preserve the above coinci- 
dence, the natural demand is to set ReII pip2 = 0. Since, 
in Eq. flAH), Ren 2 ip2 = (Ren piP2 ) 2 - (Imn piP2 ) 2 , then, 
to be precise, some mass shift survives which is equal to 



Am, 



m Pi T2 Pl ^( m2 Pl ) ( 9 P 2^ 



2(m 2 



SpiWTT 



However, even in the worse case g P2 T r7T = 0.8 (see Table HI 
where the magnitudes of the coupling constants extracted 



from the specific fits are given), this shift is estimated 
at the level of 0.1 MeV. This estimate falls within the 
errors of m pi , quoted in Table |TJ Having in mind the 
three-resonance case, we set Ren piP3 = 0. Such a type 
of justification is not applicable for the poorly studied 
resonances p2 = p(lA50) and p3 — p(1700); hence, the 
parameter 023 fixing ReiI P2P3 remains free. 

The generalization to the case of three (and any num- 
ber of) resonances pi,p2, and p^ is straightforward. The 
matrix of inverse propagators is given by Eq. (|3.4p . The 
matrix of propagators is 



G- x = 



1 



511 512 513 
.912 522 523 
513 523 533 



where 

5ii 
522 

533 

512 

513 
523 



D D — Tr 

LJ P2 ly P3 11 P2P3' 



D P1 D P2 - 

Dp 3 n plP2 
Dp 2 n plP3 

^pi^p2Pa 



2 

P1P3 

2 

P1P3 

n 



P1P3IIP2P3 1 



np 1 p 2 n P2P3 , 
npip 2 n plP3 , 



A = 



detG — Dp 1 Dp 2 Dp 3 2Hp 1 p 2 xxp 1 p 3 xxp 2 p 3 



) 11-/1, rt o II/- 



Pl P2P3 ^P2 Ai pip 3 P3 ±± pip 2 - 

(A5) 

Note that, deep in the spacelike domain, the quantity 
1/A and, as a consequence, the pion form factor have 
a pole at =87, 82, 97, and 95 GeV, when evalu- 
ated with the resonance parameters obtained from the 
fit of, respectively, SND [13], CMD-2 QJ], KLOE [3, 
and BaBaR [l5[ data. This pole is the analog of the 
famous Landau pole. 

In addition to the strong mixing between the isovector 
resonances, one should include also the isovector-isoscalar 
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Pi— w(782) mixing arising due to small G-parity breaking. 
Then, the matrix of inverse propagators can be written 
in the form 





( D Pl 


npip 2 


"■pips 






Gtot — 


— npip 3 


D P , 

— n p2 p 3 


~n P2P3 

D P3 


Hp 2 o; 
— n P3W 


(A6) 






^-pihj 




dJ J 





In this case, the pion form factor is written as follows: 

/ <?pi7T7r 



Ftt{s) — (9*1 pi 1 9*i pi 1 9~i P3 ' 9yi*>)(-* 



9p2KTT 

9p3.'t-n 
\ 9uj7nr 



(A7) 



The coupling constant g^nir describes the direct uj — > 
Tt + Tr~ transition arising due to the violation of G-parity 
conservation side-by-side with the mixing mechanism. 
However, it is known [28| that, since the 7t + tt~ channel 
dominates the p\ decay width, g UT „ is almost canceled 
in the effective u — > n + ir~ transition amplitude due to 
the compensation among imaginary parts of H PlLJ and 
the inverse p\ propagator. Indeed, allowing for both the 
mixing and direct transition, one can write the effective 
ujtttt coupling constant in the form 



o (cff) sa a 



(ReIIp lW + Tmn piW ).g Pl 

7T7I 

m i - m pi - ^ 



(ReII PlW + ilmR PlU] )g Pl im 



- m 2 pi - iy/s(T u - r pi7r7r ) 



{gu 



r 2 2 
[m u -m pi - 



ReU plLJ + i ^lmU plul 



Vs r 



pl7TTC 



m l - m Pl - iVs(Fu 



pl7T7T 



(A8) 



where ImII PlL j differs from ImII PlL j by the absence of the 
term oc g W7r7r . Hence, one can safely neglect the coupling 
constant g W7r7r . This circumstance was not properly ac- 
counted for in our earlier work, Ref. Q- The isovector- 
isoscalar type of weak mixing is essential only for the 
p(770) — w(782) system because it is enhanced due to the 
small mass difference of these resonances. See Eq. (|A8|) . 
As for other isovector-isoscalar mixings p(1450) — w(782) 
and p(1700) — w(782), there is no enhancement, due to 
the mass proximity, and one can neglect H P2aUJ . Tak- 
ing the latter assumption into account and allowing for 
the p\u mixing to first order, one can approximate the 
propagator matrix G" 1 in Eq. (|A7|) by the expression 



where the g^ and A are given by Eq. (|A5|) . The final 
approximate expression for the pion form factor F n = 
Fir(s) given by Eq. (13. ip is obtained by inserting this 
approximate expression to Eq. (|A7j) and by neglecting 
the coupling constant of the direct decay g U mr- 
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